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Abstract 


The study of the relative phase between strong and electromagnetic amplitudes 
is of great importance for understanding the dynamics of charmonium decays. The 
information of the phase can be obtained model-independently by fitting the scan 
data of some special decay channels, one of which is Y’ —^ KK. To find out 
the optimal data taking strategy for a scan experiment in the phase study of y > 
KT K" decay, the minimization process is analyzed in detail from a theoretical 
point of view. The result indicates that for one parameter fit, only one data taking 
point in the vicinity of resonance peak is sufficient to acquire the optimal precision. 
Numerical results are obtained by simulated scan fitting. Besides the results related 
to the relative phase between strong and electromagnetic amplitudes, the method 
is extended to analyze the fits of other resonance parameters, such as the mass and 
the total decay width of v/. 


1 Introduction 


The charmonium hadronic decay is mainly through two processes: the strong and 
the electromagnetic interactions. The relative phase between the strong and the electro- 


magnetic decay amplitudes is an important parameter in understanding decay dynamics. 
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Studies have been carried out for many J/q two-body decay modes: Vector-Pseudoscalar 
(VP) [1, 2], Pseudoscalar-Pseudoscalar (PP) [3, 4, 5], Vector- Vector (VV) [5] and Nucleon- 
antiNucleon (N.N) [6]. These analyzes reveal that there exists a relative orthogonal phase 
between the strong and the electromagnetic amplitudes in J/1» decays [1, 2, 3, 4, 5, 6, 7]. 
As to (W. there is also a theoretical argument which favors the +90° phase [8]. Experi- 
mentally, some analyses [9, 10, 11] based on limited VP and PP data indicate that such 
a phase is compatible with the data. Moreover, some efforts have been made to extend 
the phase study to v" decay phenomenologically [12, 13] and experimentally [14]. 

The phase study can provide valuable clue for exploring the relation between the strong 
and the electromagnetic interactions. Now with the upgraded accelerator BEPCII [15] 
and detector BESIII [16], a luminosity of 6.5 x 10??cm ?s ^! had achieved, which is the 
highest luminosity in 7-charm energy region ever existed. The 226 M J/w events, 106 M 
w’ events, and 2.9 fb ^! 4(3770) data have been collected [17], even more colossal data are 
to be collected in the forthcoming years, which gives a great opportunity to determine 
the phase between the strong and the electromagnetic amplitudes with unprecedented 
statistical precision. 

However, examining the existing determination of the relative phase, since the data are 
merely taken at one or two energy points, we find most of studies are model-dependent. 
A typical model assumption is the SU(3) symmetry in charmonium decays which supply 
additional constraint on the electromagnetic decay amplitudes in charmonium decays into 
similar final states such as VP, PP and so on. Now with high luminosity accelerator, it 
is possible to measure the phase model-independently. Since a series of energy points are 
taken for scan experiment, the model-dependent constraint no long exist. Nevertheless, 
as far as the present luminosity is concerned, the scan of charmonium exclusive decay is 
still fairly difficult due to low statistics at each energy point, especially when far from 
the resonant peak. Therefore, the optimization study for the data taking strategy is of 
great importance in order to obtain the most precise results with the limited luminosity 
(equivalently within the limited data taking time). 

Without losing generality, we focus on the mode of v/ decays to K* K- final state. 
Because, as will be shown in the next section, this decay mode can accommodate a 
comparatively simple parametrization form which is of great benefit to extract the relative 
phase. 

As far as the optimization of data taking strategy is concerned, sampling simulation 
technique was adopted for the strategy determination of 7 mass measurement [18, 19]. An 
interesting conclusion from the study is that for one parameter fit, only one data taking 
point is enough to acquire the best precision. The breakthrough of this monograph lies 
in that the minimization process is analyzed in detail from a theoretical point of view, 


which leads to the same conclusion as that of r mass measurement. As the confirmation, 


numerical results are obtained by simulated y scan fitting. Moreover, this method are 
utilized for analyzing the fits of other resonance parameters, such as the mass and the 
total decay width of y”. 


2 Theoretical Framework 


For y’ decays to Pseudoscalar-Pseudoscalar final state, the 7^7 channel is through 
electromagnetic process, the KK? through SU(3) breaking strong process, and the 
K+K- through both. Therefore, the y’ — K*K- decay is the only process which 
can be used to study the phase between strong and electromagnetic interactions by a 
scan experiment. Taken into account the continuum process, the decay amplitudes of this 
channels are parametrized as [9, 10, 20]: 


3 
Ax+k- -5,+8+ 2m, (1) 


where E. denotes the amplitude of the continuum contribution, E the electromagnet- 
ic amplitude, and {LBM the SU(3) breaking strong amplitude. They can be expressed 
explicitly as 
E. x 3 E x Hy Syr ox Ceit} p(s), (2) 
S S 2 S 
where the real parameters ¢ and C are the relative phase and the relative strength between 


the strong and the electromagnetic amplitudes, and B(s) is defined as [9] 


B() = veel (3) 
s— M? +iMT: 
where ys is the center of mass energy, o is the QED fine structure constant, M and T; 
are the mass and total width of Y’, Fee is the partial width to e*e~ final state. 
'The Born order cross section for this channel reads 


Ara? 
o Born = SE [1 + 22:(C4B(s)) + |CGB(s)?] x [Free x- (5)? Pe - (8), (4) 
where Cy = 14-Ce'?; Fg+g-(s) is the form factor, which is usually written as Pre (s) = 


fr+x-/s with fx+x- being a constant; Px+x-(s) = oak is the phase space factor, with 


E 
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where gx is the momentum of K* or K^, my is the nominal mass of K meson. 
In actual experiment, the effect of Initial State Radiation (ISR) is considered through 
an integral [21] 


Wm i, " AE meee (5) 


where F(x, s) is the structure function which can be calculated to an accuracy of 0.1% [21, 
22, 23]. 

In addition, another important experimental effect, the energy spread of et and e^ 
must also be taken into consideration. Finally, the experimentally observed cross section 


is expressed as |24, 25] 


sss) = | "d gas sts atus. uo), (6) 


0 
where G(s’, /s) is a Gaussian distribution 
l _(Ws-va)? 
€ 282 — 
V27A 


Here A indicates the energy spread of the collision beams. 


QLY S, V/s) = 
Some parameter values for the numerical calculation in the following sections are 
articulated in Table 1. 


Table 1: Some parameter values for numerical calculation. The quantity with x will be 


set as a free fitting parameter for the corresponding study. 


Quantity numerical value Remark 


xM 3.68609 GeV Ref. [26 
xT: 304 keV Ref. [26 
Lom 2.35 keV Ref. |26 
Mk 493.677 MeV Ref. [26 
P e 0.9 GeV? Ref. [27 
A 1.3 MeV Ref. [28 
xo 90? Ref. [29 
C 2.5 Ref. [29 


3 Minimization Analysis 


For a scan experiment, several points, say totally N,, points, need to be taken in the 
vicinity of the resonance (in this monograph the resonance is Y’). The estimator is usually 


constructed as a X? function [30]: 


2 (Ne = Laes: 
X E ES (ANots)2 H (7) 
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where N? and AN?* = ,/N°s are the observed number of events and its error at 
the i-th point, L; the corresponding luminosity, e; the selection efficiency, and o; the 
theoretical cross section that is desp in Eq. (6). The fitting parameters (relative phase, 
strength, etc.) are contained in c;, and these parameters and the corresponding errors 
can be extracted by minimizing the x? function defined in Eq. (7). In the following 
analyses, only concerned is one free fitting parameter, the relative phase between strong 
and electromagnetic amplitudes, that is, 0. 


bs 


If we denote the observed cross section measured at energy point i as 02^, and rewrite 


Me = Lo S (8) 


Eq. (7) can be recast as 


j= Y g^ — gl YA Aem — 0) (9) 


Here x? is replaced with f for simplicity, and the following relations are utilized : 
Li = xilo , y» =l; (10) 


where Lo is the total luminosity (corresponding to the finite total data taking time) and 
x; is the fraction of luminosity at the i-th energy point. Moreover, Z: is supposed to be 
the same at all points (e = 50% is used for numerical calculation), which is a fairly good 
approximation for the scan of narrow resonances, such as J/w and y”. 

In the light of Eq. (9), the first and second order derivatives of the function f to @ 


can be derived as 


af S s usa 00% 
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The experimentally concerned cross section functions are generally smooth enough, 


which can be approximated by polynomial functions. Therefore, the first and second 
order derivatives of these functions are also smooth enough (refer to Fig. 1). Under such 
case, we argue that the second term in Eq. (12) could be neglected. When the fitting 
process finishes, c; in Eq. (11) and Eq. (12) could be considered as the true value of the 


obs 


cross section at energy point 7. As we assumed previously, 92^ is the experimentally 


obs 


measured cross section at energy point i, then (02^ — o;) could be considered as a random 


variable which satisfy a Gaussian distribution with mean as 0 and deviation as Ag? (the 


obs ). 


error of af As a conservative estimation, we assume the relative error of cross section 


0.0.07 / pb 


Figure 1: (a) Comparison between the cross section and its first and second order deriva- 
tive. The solid line shows the cross section (c). The dashed line shows the first order 
derivative (o^). The dotted-dashed line shows the second order derivative (o"). (b) Com- 
parison between first and second term inside the bracket in Eq. (12). The solid line only 
shows the value of the first term (T1). The dotted line contains contributions both from 
the first term and the second term (T; + Tə). 


measurement is 10%, which means Ag% = 0.1 - he 


The expectation of the second 
term inside the sum in Eq. (12) could be calculated by using sampling method. The 
comparison of these two terms are shown in Fig. 1, from which the second term is quite 
small compared with the first one, therefore its effect can be neglected. Now Eq. (12) 


becomes 


Oe f (13) 
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zi (ci 
= 2 LoE » = == 
-— — 0; V 0¢ 
$—Ó9 i t g=¢* 
where ¢* is the fitting result of the relative phase and o7 is the theoretical cross section 


at d = ¢”*. 
The fitting error of à can be evaluated as [30] 


(14) 


According to Eq. (14), the maximum of the second order derivative of fitting function 


yields the minimum of fitting error. Define a new function g as 


1 £3 
* a 
7i 6 o=¢* 
where the subscript i denotes the value of g at the i-th energy point. Then Eq. (13) 
becomes 
az =2Le x (16) 
0g? m — 0 iGi- 
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Notice that $7; x; = 1, it is readily to obtain the following inequalities 
Ümin = om 32) Gina < X rigi = DD Li) Jmaz = Imaz, (17) 


where gmin (Gmar) is the minimum (maximum) value of g within the energy region con- 
cerned. To get maximum 5 7, z;g;, only one data taking point is sufficient and it should 


be located at the energy point which renders g maximum. 


4 Numerical results 


To reinforce the preceding conclusion, the simulated scan data are fit to get numerical 
results. In this procedure, great many times of fitting need to be performed, where the 
large number of calculations must be carried out for the observed cross section. Unfor- 
tunately, two nested integrations of the observed cross section, which take into account 
both the ISR correction and beam energy spread effect, take so much time that any actual 
optimization fitting becomes impractical. In a recent study [31], using the generalized lin- 
ear regression approach, a complex energy-dependent factor is approximated by a linear 
function of energy. Taking advantage of this simplification, the integration of ISR correc- 
tion can be performed and an analytical expression with accuracy at the level of 1% is 
obtained. Then, the original two-fold integral is simplified into a one-fold integral, which 
reduces the total computing time by two orders of magnitude. In the following studies, 


the simplified observed cross section formulas are adopted to acquire all numerical results. 


4.1 Relative Phase 


Considering the parameter to be analyzed is the relative phase between strong and 
electromagnetic amplitude of y’ decay, the distribution of gẹ and the fitting error Ej on 
energy region when $* = 90° is shown in Fig. 2, according to which, at the energy point 
3.6868 GeV, the value of function gg reaches its maximum while E, reaches its minimum. 
In the vicinity of 3.684 GeV, the gg value is very small and the corresponding Ey is quite 
large. So this point (3.684 GeV) should be avoided in the scan experiment!. 

By fixing the energy point to 3.6868 GeV, the error obtained from fitting and computed 
by Eq. (14) is shown in Fig. 3. Just as expected, the higher the luminosity, the smaller 
the error. Moreover, the fitting and computing values of error are so consist with each 
other that it is hardly to distinguish them in Fig. 3 (a). To exhibit the details, the relative 
difference of Eg, that is, 

E,(computing) — E;(fitting) 


[e 
i E4(computing) 


'To validate this result, the sampling technique is used to check the data taking point distribution 
and the fitting error. Details about the sampling technique can be found in Appendix. 
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Figure 2: The distributions of gẹ and E; for ¢* = 90°. 


is shown in Fig. 3 (b). 

If the relative phase @ variates, the optimal position of energy will change correspond- 
ingly. Table 2 lists the optimal values of energy position for some special phase angles. 
According to these information, the values of the optimal energy points arrange from 
3.686 GeV ~ 3.687 GeV, nearly within the scope of 1 MeV. 


Table 2: The optimal data taking position for different relative phase @ 
OL") 0 45 90 135 180 270 
Optimal point (GeV) | 3.68604 | 3.68700 | 3.68680 | 3.68706 | 3.68648 | 3.68672 


4.2 Other parameters 


The minimization analysis discussed in Section 3 is applicable to any parameter we 
are concerned with. As long as the variable ¢ is replaced with the parameter to be 
analyzed, all aforementioned deductions are valid. In study that followed, we perform the 


optimization for the two interested resonance parameters, mass and total width. 


4.2.1 7’ Mass 


The similar analyses are performed for the mass of wv" resonance, and results are 


displayed in Fig. 6, where four curves corresponding to different relative phases, ó — 


Computing Error Relative Difference 


Fitting Error --------~ 


E, / rad 


0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 “0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
Luminosity / pb! Luminosity / pb?! 


a 


Figure 3: (a) Comparison of the equivalent error for different luminosity by fixing the 
energy point at 3.6868 GeV; (b) Relative difference Ry at 3.6868 GeV for different lumi- 


nosity. 


0°, 90°, 180°, and 270°. There are two new features for the optimization of mass parameter. 
Firstly, the energy position for smallest Ej is at 3.6845 GeV for ó = 0°, 90°, and 270°; 
but at 3.6874 GeV for o = 180°. Secondly, two energy positions should be avoided due to 
larger values of Em. One is around 3.686 GeV, the other is near 3.68 GeV for ¢ = 0?,90?, 
and 270°; while near 3.69 GeV for ¢ = 180°. 

To understand the heterogeneous behavior of curve for ¢ = 180? from the other ones, 
we take the curve for ó = 90° as a representative, and show the first order derivatives 


of cross sections in Fig. 5, where the left and right rows correspond to 90° and 180°, 


R "n 


respectively. The total cross section ezp is divided into three parts, ie. ohp Texp, and 


012, as we did in Ref. [31]. The corresponding derivatives of them are shown sequently 


in Fig. 5. Investigation of those figures indicates that the crucial role for the different 


I1 
exp? 


which is opposite to each other. Furthermore, if we scrutinize the equations in Section 
3 of Ref. [31], the sign of the derivative of o/1, is determined by the coefficient A> = 


exp 


6(Tee/a) - (1 + C cose), which is positive when ¢ = 90° and negative when à = 180° 
for C = 2.5. The switching point is at ¢ = cos !(—1/C) ~ 113.5°. 


behavior between 90°- and 180°-curves is played by the o2,.,, the derivative variation of 


4.2.2 wy’ Total Width 


The optimization results for the total width of v" resonance are shown in Fig. 6, where 
three curves corresponding to distinctive beam energy spreads, A = 1.1, 1.3, 1.5 MeV, are 
presented. With the enhancement of A, the position of minimum error, Erp, shifts a 
little bit rightward along the abscissa. For all circumstances, the energy position for the 


maximum Ep is near 3.68 GeV while for the minimum error is around 3.686 GeV. 
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Figure 4: The results of analyzing v/' mass for relative phase ¢ = 0°, 90°, 180? and 270°. 
(a) is the gm value and (b) is the fitting error Em. 


5 Conclusion and Discussion 


In this paper, one-parameter-and-one-point conclusion is demonstrated through a the- 
oretical analysis of minimization process instead of the sampling simulation as we did 
before. As far as the phase study is concerned, for the y’ — K+K- process, the optimal 
data taking point is determined to at 3.6868 GeV which is near the peak of 17 nominal 
mass. The same method is also used to acquire the optimal point for other resonance 
parameters, such as the mass and the total width of y”. 

In principle, the idea put forth in Section 3 could be extended for multi-parameter 
optimization. Formally speaking, the vector and matrix quantities, would be adopted for 
the corresponding analysis, say, the second order derivative of one parameter are to be 
replaced by Hessian matrix, a matrix of second order derivative for all parameters. 

However, there are some problems not easily to be settled. The most prominent one 
is how to define “optimal”. In one parameter scenario, the optimal data taking point is 
the one which could make the fitting error of the parameter reaches its minimum. But 
for multi-parameters, there are many options: the sum of relative fitting errors of all 
parameters reaches its minimum; the merely fitting error of one major parameter reaches 
its minimum while others do not. Different options lead to distinctive results. All this 


makes the situation more complicated and is left to the study in the further. 


Appendix: Sampling Technique Methodology 


Suppose there are Np data taking points in experiment, and the theoretical number 


of events in i-th energy point could be calculated as 


Nite — Lj-0;- 6,42 1,2,..., Ng, (18) 
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Figure 5: The derivative distributions for different parts of the cross section. (a) — (d) is 
for ¢ = 90° and (e) - (h) is for ¢ = 180°. (a)(e): ee. (b)(f): o£, (c)(g): ell. (d)(h): 


exp’ exp’ 
I2 


mn» 


11 


A=1.1MeV 


gp /pb Gev'! 
Ep/ GeV 


Figure 6: The results of analyzing Y" total width for A = 1.1 MeV, 1.3 MeV and 1.5 
MeV. (a) is the gr value and (b) is the fitting error Er. 
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Figure 7: (a) Scan error distribution for five consecutive points; (b) The relative difference 


between five consecutive points and one point scan. 


where L; is the integrated luminosity in i-th energy point, e; is the event selection effi- 
ciency. 
In sampling technique, the experimentally observed number of events and its error 
could be taken as 
AN? = (Nthey1/2 (19) 


Ne = Nee + £ g AN, (20) 


where £ is a random number which satisfy Gaussian distribution. 

Using the observed event number and its error calculated above, the parameter we 
interested in (in this paper, the relative phase) could be fitted by finding the minimum of 
Eq. (7). 

By repeating this process, a large number of observed event number and error could 
be generated and so does the fitting parameter and its error. We can compare these errors 
with the results obtained from the method we just developed. 

In section 4.1, we scan through the energy region using one data taking point. To 


check this result, five consecutive energy points are used in energy scan. The fitting 


12 


error of five consecutive points scan versus the central energy point is shown in Fig. 7. 


This result is similar with one point scan. Their difference is also shown in Fig. 7. The 


difference between the two scan schemes is generally at the level of one per mille except 


for the points around 3.684 GeV, where the variation of error curve is rather rapidly. As 


to the five-point scheme, there is at least one point within the region with comparatively 


large error. 
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